Oscillations of atomic fermions in a one dimensional optical lattice 
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A semiclassical model is used to investigate oscillations of atomic fermions in a combined magnetic 
trap and one dimensional optical lattice potential following axial displacement of the trap. The 
oscillations are shown to have a characteristic small amplitude, damped behavior in the coUisionless 
regime. The presence of a separatrix in the semiclassical Brillouin zone phase space is predicted and 
shown to produce a strongly asymmetric phase space distribution function. 
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I. INTRODUCTION 

The recent advances in the cooling of atomic Bose and 
Fermi gases to degeneracy, and the simplicity with which 
the gas can be confined in a periodic optical lattice poten- 
tial provide a new point of departure for quantum trans- 
ort studies of atomic bosons Q, H, H, 1^ P and fermions 



In particular, neutral atomic fermions form a 
nearly ideal gas at nanokelvin temperatures, and by con- 
trast to conduction electrons in solids, can traverse the 
dimensions of the Brillouin zone without collision. 

In a recent experiment a gas of fermionic potassium 
40 atoms was evaporatively cooled in a combined mag- 
netic and optical lattice potential, leading to a degener- 
ate cloud at a temperature of a fraction, approximately, 
one third the Fermi temperature Tp for the system of N 
atoms ^ . Specifically, the experiment involved a config- 
uration with a one dimensional axial lattice potential pro- 
duced by a laser standing wave field, with tight magnetic 
confinement in the transverse plane, but a weaker mag- 
netic trap in the direction coaxial with the lattice thereby 
emphasizing the latter's role. By displacing the center of 
the axial magnetic field with respect to the lattice, dipo- 
lar oscillations of the Fermi cloud in the lattice poten- 
tial were investigated. The principal observations were 
that the amplitude of the oscillations were significantly 
smaller than the trap displacement, and were damped 
over several periods of the motion. 

In this paper we show that these observations can be 
explained by a simple semiclassical theory of coUisionless 
atomic fermion motion in the combined magneto-optical 
potential. This confirms the qualitative physical picture 
discussed in ref.0. The effects of fermion statistics are 
felt only in the initial conditions. A broad distribution of 
occupied quasi momentum states with large radial quan- 
tum numbers is created by evaporative cooling on ac- 
count of the Pauli exclusion principle. In a semiclassical 
picture the lattice transport dynamics can in principle 
be determined by following the particle trajectories in a 
Brillouin zone phase space consisting of atomic position 
and quasi momentum. In the problem of interest here, 
we are able to do this explicitly in the reduced two di- 
mensional phase space of axial components of position 
and quasi momentum. The effects of radial motion are 



readily accounted for by an appropriate summation over 
radial quantum numbers. 

The remainder of this paper is organized as follows. 
In section II we develop a theoretical model for the dy- 
namics of the degenerate cloud in the confined magnetic 
and optical lattice potentials. In section III we discuss 
the semiclassical distribution function on the axial phase 
space, and illustrate how the trap displacement influences 
the evolution of the distribution function, and the center 
of mass motion of the cloud. In section IV we summarize 
our conclusions. 



II. SEMICLASSICAL ATOMIC PHASE SPACE 
TRAJECTORIES 

The Hamiltonian for single particle motion in a com- 
bined optical lattice and magnetic confining potential 
may be usefully written in the following form 



H = Ho + Hi + V{z), 



(1) 



where Hq = {pi + j3^)/(2m) + muj'^{x'^ + y'^)/2, Hi = 
pl/2m+ Vosin'^kLZ and V{z) = The Hamil- 

tonian Hq governs the radial harmonic motion produced 
by the magnetic trap, while V{z) is the residual part 
of the magnetic confinement in the axial direction. For 
tight radial confinement >> uja- The Hamiltonian 
Hi{z) = Hi(z + a) describes one dimensional motion in 
an optical lattice with period a = Ai/2, where Xl is the 
optical wavelength and kL = 27r/AL is the wavenumber. 
The lattice depth Vq = sEr is conventionally parameter- 
ized by a dimensionless multiple s, of the atomic recoil 
energy Er — fi?k\/ {2m). 

We are here interested in the following specific ex- 
perimental situation. A single component Fermi gas is 
cooled to degeneracy in the combined magnetic and op- 
tical lattice potential to well below the Fermi tempera- 
ture Tp. The axial potential is then rapidly displaced 
V{z) V{z — d), and the degenerate gas performs oscil- 
lations in the lattice. These oscillations were observed by 
waiting for some time before switching off both magnetic 
and trap potentials and imaging the cloud after a period 
of ballistic expansion |^ . 
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The one dimensional lattice Hamiltonian Hi has Bloch 
eigenfunctions ipaq{z) and eigenvalues Sail) = £ai<l + K) 
with K — (27r/a) x integer, any reciprocal lattice vec- 
tor. Here hq is the axial quasi momentum associated 
with translations z — > z + a in the periodic lattice, and 
a is the integer axial band index. We employ a semi- 
classical approach to treat the transport dynamics of the 
atomic fermions in the axial direction. The method is 
well suited to the description of the coUisionless regime 
relevant to single component ultracold fermions where 
s wave scattering is prohibited by the antisymmetry of 
the two-body wave function and p-wave scattering is 
strongly suppressed. The energy bands implicitly take 
into account the axial lattice potential responsible for 
large changes in the physical momentum of the moving 
particle, but which conserve the quasi momentum hq. We 
assume, consistent with ref. , that only the fundamen- 
tal axial band a = 1 is populated at the temperatures 
of interest, and therefore we will suppress the axial band 
index. In order to correctly describe quantum mechani- 
cal motion in the plane transverse to the lattice, radial 
quantum numbers will arise later in the theoretical de- 
velopment. 

In the semiclassical theory, motion of an atom is con- 
sidered to be restricted to a given energy band, and only 
external forces other than the lattice, in this case due to 
the magnetic trap, can change its quasi momentum. It 
follows that the axial motion of a fiducial atom, in re- 
sponse to the axial displacement d of the trap, is found 
by solving the equations of motion 



hq{t) 

m 



-mujl{z{t) - d) 
1 d£{q) 
h dq 



(2) 



We note that these equations are independent of any 
radial quantum numbers associated with motion in the 
transverse xy plane. 

As a simple example consider motion confined to 
near the bottom of the fundamental axial band, but 
with arbitrary radial excitation. We may write £{q) — 
S^(7^/(2me//) -I- constant, with me// the associated effec- 
tive mass. The equations of motion reduce to 



d'^zjt) 
dfi 



-if {z{t)~d), 



(3) 



describing harmonic oscillations about the displaced po- 
sition d with renormalized frequency = ^m/meff uja 
\l(f\. The harmonic approximation predicts large ampli- 
tude dipole oscillations but cannot correctly describe the 
orbits of occupied fermionic states with relatively large 
axial quasi momentum far from the band minimum. On 
the other hand a Bose condensate would have a narrow 
quasi momentum distribution Q. The simple harmonic 
theory indicates that one should expect large amplitude 
dipole oscillations in this has been observed in 

experiments 

An ideal gas model with a narrow 
quasi momentum distribution is, however, not sufficient 



to correctly describe superfluid effects associated with an 
interacting Bose condensate (12j . 

In order to treat non harmonic phase space orbits we 
note that any axial energy band must satisfy the peri- 
odicity requirement £{q) = £{q + K). For example, the 
axial energy band 



irq 

2k~L 



(4) 



behaves like li^q^ / {2m^f f) for small g, but also satisfies 
the periodicity requirement with period 2fci. Thus all 
physically distinct states of a single fermionic atom lie in 
the range — /cl < q < kL- This energy band has the ad- 
ditional advantage that it produces pendulum dynamics 
and the model can be solved analytically. In this case 
the equations of motion become 



hq{t) = —muj'^{z{t) — d) 

m = 



hk^ f TTq(t) 
sm ' 



(5) 
(6) 



Eliminating z and defining the dimensionless quasi mo- 
mentum Q = nq/ (2A:i) we get the conservation equation 



d_ 

dt 



0. 



(7) 



The first integral may be written Q{t)^ + il^sin^Q{t) = 
ri^M where M is fixed for a given orbit. The solution is 
given by 



sinQit) = VMsn (ifii + (j)M\M) (8) 
where the initial phase 



<I,M = u{sinQ{0)/VM\M) 



(9) 



sinQ{t) = sn (±VAfl7i -|- 0i/jv/|l/Af 



is given in terms of the Legendre elliptic func- 
tion of the first kind u(a;|m). Using the identities 
sn{x\M) = sn(VMa;|l/M)/\/M and \/Mu(sm(/)|A/) = 
u{-s/M sin(f)\\ / M) ^ we can also write 

(10) 

where the initial phase 4>i/m — 't'M^M is given by 

4>i/M^u{sinQmi/M). (11) 

We present these results in terms of elliptic functions in 
which the amplitude argument, either M or 1/M, is in 
the range [0, 1] . This is useful for the comparison of orbits 
with M < 1 and 1/AI < 1. As we discuss further below, 
these orbits are analogous to the oscillating and rotating 
orbits, respectively, of a mechanical pendulum. 

Differentiation of the result for sinQit) with respect to 
time yields 

Q{t) = ±Vt\fMcn{±Q.t + (t)M\M), (12) 
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which in turn gives the atomic trajectories 

z{t)-d = T^VMcn{±nt + (f)M\M). (13) 
or, using cn{x\A4) = dn{VMx\l/M) 



z{t) -d^ TiVMdn{±VMnt + 0i/m|1/M). 

The value of M depends on the trap displacement and 
satisfies the "energy" conservation condition, 



(z(0) - df + esin^Q{0) = fM. 
where we define the length scale 



eff \mU!a 



(14) 



(15) 



such that £{q) = eosin^Q with eo — muj'^£'^/2. For 
M < 1, the orbits in the z — q phase space are closed, 
whereas for M > 1 they are open. The separatrix or- 
bit M = 1 has homoclinic fixed points on the Brillouin 
zone boundary, analogous to the infinite period motion 
of a physical pendulum between vertically upwards ini- 
tial and final positions, sinQ{t) = tanh(±nt + 0i) and 
{z{t) — d)/£ = ^sech{±rtt + (pi). The qualitatively dif- 
ferent character of these orbits is responsible for a non- 
trivial dependence of the dipole oscillations on the initial 
trap displacement. When the trap is displaced z ^ z — d, 
the peak of the phase space distribution function remains 
at z = 0, but the separatrix is given by Ea. (|14|) with 
M — 1, and is centered in the region z > 0. Subse- 
quently the atoms lying inside the separatrix will perform 
closed orbits, while those left outside remain outside, and 
perform small amplitude spatial oscillations as they tra- 
verse the Brillouin zone. The relative importance of the 
two sets of orbits in the overall motion of the cloud is 
determined by the fraction of atoms initially inside and 
outside the separatrix. In turn, this is governed by the 
magnitude of the initial displacement d relative to i. As 
we will illustrate in the following section, the axial distri- 
bution function can develop strong asymmetry for large 
trap displacements, due to the separatrix structure of the 
Brillouin zone phase space. 



III. DYNAMICS OF THE AXIAL ATOMIC 
DISTRIBUTION FUNCTION 

The initial thermal distribution of atomic fermions pro- 
duced by evaporative cooling in the combined magnetic 
trap and optical lattice potentials, can be written in 
terms of a semiclassical distribution function for motion 
in the axial direction with the radial harmonic oscillator 
quantum numbers and Uy. Since the radial poten- 
tial is assumed to be azimuthally symmetric, only the 
combination + riy arises. Recall that we assume the 
temperature is sufficiently low that only the lowest axial 
band is populated 9|. We are therefore left to consider 



a set of axial distribution functions with different radial 
quantum numbers. Specifically, the initial axial distri- 
bution function is given by f^°\z, q) = X]^=o fn}{z, q), 
where 



/(°)(z,9) = (n„ + l) 



^l3(n^huj^+£(q) + V(z)-tJ.) 



1 



(16) 

with (3 — l/{kBT), rir = 0,1,2,... the radial quantum 
number and fj, the chemical potential. Note that, due 
to the lattice, this is a function of z and quasi momen- 
tum hq, as opposed to physical momentum. The factor 
{rir + 1) is the degeneracy of the radial band with energy 
fiWrTLr abovc the ground state rir = 0. The distribution 
function is assumed not to change if the axial trap is dis- 
placed rapidly. The subsequent coUisionless motion of 
the cloud causes the distribution to change according to 
the Boltzmann equation 



dfn^ , d£{q) dfn,. 



dt 



dhq dz 



with solution 



muu^ {z — d) 



/„^(z,9,<)-/(°)(z(t),g(t)) 



dfn 



0, (17) 



(18) 



Using the previous analysis above we have that, for M < 
1, 

sinQ{0)cnntdnnt - ^^^"j'"'^ cosQ(0)dn»t 

_ (z(0) - d)cnnt + lsinQ{0)cosQ{Q)snVltdnQ.t 
1 - sin^Q{Q)sv?nt 



z{t) 



where we note that (in(0jv/|A/) = cosQ(O) > in the 
first Brillouin zone, and snU,t = sn{rit\M); similarly for 
cn and dn. For A/ > 1, the following expressions hold 



sinQit) 



z{t) - d 



(20) 



sinQ{0)cnujtdnujt — ^^^^=^cosQ{0)snujt 

1 - sin^Q{0)sn'^ujt/M 
(z(0) — d)dnujt + -^^sinQ{0)cosQ{0)snujtcnujt 
1 - sin^Q{0)sn^ujt/M 



where we note cn{4'i/]^j\l/M) = cosQ{0) and sncut — 
sn{VMnt\l/M), etc. 

The axial number density of fermions in radial band 
Tir at time t, is given by integrating over the quasi mo- 
mentum of a single Brillouin zone 

f^^ dq 
J-kL 27r 

and the total axial column density n(z, t) = 
^^^Q ri„^(z, t). Similarly, the distribution of quasi mo- 
mentum can be computed by summation of the distri- 
bution functions for the individual radial energy bands, 
rinAl^t) = I dzf„^{z,q,t). 

In Fig.l, we show the initial spatial density as a func- 
tion of radial quantum number and z at a temperature 



FIG. 1: The initial atomic density is shown as a function of 
the radial quantum number n,.. The total axial density n(z, 0) 
is found by summation over rir- Parameters are given in the 
body of the text. 
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FIG. 2: The initial axial atomic distribution function in the 
z-Q phase space for d = 15^m. The oval shaped contour is 
the separatrix across which particles may not cross during the 
coUisionless dynamics. Other parameters are the same as in 
Fig. El 



T « O.STp. We have taken the experimental parameters 
from ref. 0, except that we use the analytical model 
energy band discussed in the last section rather than 
the numerical result appropriate to an optical standing 
wave. Hence our results are not intended for quantita- 
tive comparison with the experiment. However, since 
the qualitative features of the model energy band are the 
same as those shown in ref. , its predictions should be 
in qualitative agreement with the experimental observa- 
tions. The parameters employed in the figures are given 
here mainly in units of temperature, therefore energies 
are recovered by multiplying by Boltzmann's constant. 
We consider potassium-40 atoms cooled in a magnetic 
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FIG. 3: The axial atomic distribution function at Qt — 15, in 
the z-Q phase space for d = 15/^m. The region of high density 
lies inside the separatrix, but has been stretched out over the 
region. Other parameters are the same as in Fig. 
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FIG. 4: Small amplitude damped oscillation of the axial com- 
ponent of the cloud center of mass near the undisplaced center 
of the cloud z = 0. The position of the displaced trap center 
is shown by the dashed line. Other parameters are the same 
as in Fig.Q 



trap with radial and axial frequencies tOr — 15.2nK and 
LUa = 1.15nK, respectively. The wavelength of the in- 
cident light is taken to be Al = 754nm. The chemi- 
cal potential fi = 430nK, and /3 = 6.98(MK)~\ For 
our model axial energy band, assuming rrieff = m for 
simplicity, this gives the length scale £ ~ 55.5/im and 
band depth eo ~ 170nK. A density plot of the initial 
distribution function is shown in Fig. 2, along with the 
separatrix trajectory in the phase space following trap 
displacement. Orbits inside and outside the separatrix 
follow its contour precisely as occurs in the case of the 
mechanical pendulum. In this example we assume a trap 
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FIG. 5: The initial atomic distribution function in the z-Q 
phase space for d — 45/im. The region of high density tra- 
verses the phase space separatrix. Other parameters are the 
same as in Fig.0 



1000, 




FIG. 6: The atomic distribution function in the z-Q phase 
space a,t Q,tf — 15 for d = Abfim. The complicated structure 
reflects the open and closed orbits in the first Brillouin zone, 
which lead to a strong asymmetry. Other parameters are the 
same as in Fig.0 



displacement d — 15/.tm much less than i, so that the 
bulk of the atomic distribution lies inside the separatrix, 
and therefore most though not all particles will undergo 
bound orbits in the first Brillouin zone. 

The dynamics of the particle distribution in each radial 
band can be calculated separately, as explained above, 
and then summed over the radial quantum number Ur 
to get the total axial distribution function. However, the 
latter summation is unnecessary for a qualitative under- 
standing of the dynamics. Each radial band distribu- 
tion function has the same dynamics, because the phase 
space trajectories {z(t),q{t)) are independent of the ra- 



dial quantum number n^. In particular the damped mo- 
tion of the cloud center of mass is evident in the distri- 
bution function /„^(z,g, i) dynamics, for a fixed radial 
quantum number n^. In Fig|31 we show the total axial 
distribution which evolves from the initial distribution 
shown in Fig[21at time ^tf = 15. During this time, the 
center of mass of each radial band density undergoes a 
similar damped dipole oscillation leading to the center of 
mass dynamics for the cloud illustrated in Fig^J We note 
that the amplitude of the oscillation is small compared 
to the displacement of the trap center, in qualitative 
agreement with the observations of ref.[3. The damp- 
ing is caused by the stretching of the distribution func- 
tion along the phase space trajectories, not by collisions, 
dissipation or the influence of different radial orbits. It 
is straightforward to check, and is indeed obvious in the 
limit of a very narrow distribution in (z, Q) phase space, 
that the damping is correspondingly reduced. However, 
a Fermi distribution must necessarily have a broad dis- 
tribution of Q states on account of the Pauli exclusion 
principle. Damping is therefore a feature of coUisionless 
Fermi transport dynamics. A narrow quasi momentum 
distribution is predicted also by the Gross-Pitaevski dy- 
namics of a superfluid Bose condensate in a lattice, and 
in this case undamped oscillations have been observed 

illi. 

The non harmonic orbits of large axial quasi momen- 
tum which are necessarily populated in a Fermi gas cause 
the stretching of the distribution function both inside 
and outside the separatrix. As noted above, each radial 
band executes a similar damped motio in phase with the 
others. In Fig|Sl and FigEl we illustrate the dynamics 
of the axial distribution function for a larger trap dis- 
placement d = 45/zm, a value which is comparable to the 
separatrix distance £ « 55/im. In this case a fraction 
approaching one half of the particles now follow open 
orbits in the first Brillouin zone. The total distribution 
function develops strong whorls inside the separatrix and 
spatial oscillations occur in the exterior. The basic left- 
right asymmetry survives the integration over the quasi 
momentum required to determine the axial spatial den- 
sity, but the center of mass motion is qualitatively similar 
to that illustrated in Fig0] Damped oscillations of the 
Fermi gas alone cannot be viewed as strong evidence of 
a separatrix structure in the Brillouin zone phase space, 
since both closed and open orbits produce small ampli- 
tude center of mass oscillations. However, the presence 
of the separatrix could well be observable experimentally 
if the amplitude of center of mass oscillations were sys- 
tematically studied as a function of trap displacement d. 
Alternatively, direct measurements of the axial distribu- 
tion function or spatial density would be most useful. 



IV. CONCLUSION 

We have shown that a simple semiclassical theory of 
the axial motion of a gas of atomic fermions, without any 
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collisional interactions, produces transport dynamics in 
qualitative agreement with recent experimental observa- 
tions The gas executes damped oscillations with an 
amplitude small compared to the axial magnetic trap dis- 
placement. Wc also point out the presence of a separatrix 
in the Brillouin zone phase space which produces strong 
asymmetries in the axial phase space distribution and 
axial density functions. It would be interesting to inves- 
tigate this structure further experimentally, for example, 
by looking for non-monotonic behavior of the transport 
properties as a function of the initial trap displacement 



d. 
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